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Abstract. We study interactions between Skolem Arithmetic and cer¬ 
tain classes of Constraint Satisfaction Problems (CSPs). We revisit re¬ 
sults of Glafier et al. in the context of CSPs and settle the major 
open question from that paper, finding a certain satisfaction problem on 
circuits to be decidable. This we prove using the decidability of Skolem 
Arithmetic. We continue by studying first-order expansions of Skolem 
Arithmetic without constants, (N; x), as CSPs. We find already here a 
rich landscape of problems with non-trivial instances that are in P as 
well as those that are NP-complete. 

1 Introduction 

A constraint satisfaction problem (CSP) is a computational problem in which the 
input consists of a finite set of variables and a finite set of constraints, and where 
the question is whether there exists a mapping from the variables to some fixed 
domain such that all the constraints are satisfied. When the domain is finite, 
and arbitrary constraints are permitted in the input, the CSP is NP-complete. 
However, when only constraints from a restricted set of relations are allowed 
in the input, it can be possible to solve the CSP in polynomial time. The set 
of relations that is allowed to formulate the constraints in the input is often 
called the constraint language. The question which constraint languages give 
rise to polynomial-time solvable CSPs has been the topic of intensive research 
over the past years. It has been conjectured by Feder and Vardi [14] that CSPs 
for constraint languages over finite domains have a complexity dichotomy: they 
are either in P or NP-complete. This conjecture remains unsettled, although 
dichotomy is now known on substantial classes (e.g. structures of size < smm 
and smooth digraphs mm)- Various methods, combinatorial (graph-theoretic), 
logical and universal-algebraic have been brought to bear on this classification 
project, with many remarkable consequences. A conjectured delineation for the 
dichotomy was given in the algebraic language in m- 
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By now the literature on infinite-domain CSPs is also beginning to mature. 
Here the complexity can be much higher (e.g. undecidable) but on natural classes 
there is often the potential for structured classifications, and this has proved to 
be the case for reducts of, e.g. the rationals with order [3], the random (Rado) 
graph [7] and the integers with successor ; as well as first-order (fo) expansions 
of linear program feasibility [3]. Skolem Arithmetic, which we take here to be the 
non-negative integers with multiplication (and possibly constants), represents a 
perfect candidate for continuation in this vein. These natural classes have the 
property that their CSPs sit in NP and a topic of recent interest for the second 
and third authors has been natural CSPs sitting in higher complexity classes. 

Meanwhile, a literature existed on satisfiability of circuit problems over sets 
of integers involving work of the first author m, itself continuing a line of 
investieation beeun in [36| and Dursued in |38I39I27|. The oroblems in m can 
be seen as variants of certain functional CSPs whose domain is all singleton sets 
of the non-negative integers and whose relations are set operations of the form: 
complement, intersection, union, addition and multiplication (the latter two are 
defined set-wise, e.g. A x B := {ab : a € A A b € B}). An open problem was 
the complexity of the problem when the permitted set operators were precisely 
complement, intersection, union and multiplication. In this paper we resolve that 
this problem is in fact decidable, indeed in triple exponential space. We prove 
this result by using the decidability of the theory of Skolem Arithmetic with 
constants. In studying this problem we are able to bring to light existing results 
of [13 as results about their related CSPs, providing natural examples with 
interesting super-NP complexities. 

In the second part of the paper, Skolem Arithmetic takes centre stage as we 
initiate the study of the computational complexity of the CSPs of its reducts, 
i.e. those constraint languages whose relations have a fo-definition in (N; x). 
CSP(N; x) is in P, indeed it is trivial. The object therefore of our early study is 
its fo-expansions. We show that CSP(N; -b, is NP-complete, as is CSP(N; x, c) 
for each c > 1. We further show that CSP(N; x,{7) is NP-complete when U is 
any non-empty set of integers greater than 1 such that each has a prime factor 
p, for some prime p, but omits the factor Clearly, CSP(N; x,U) is in P 
(and is trivial) if U contains 0 or 1. As a counterpoint to our NP-hardness 
results, we prove that CSP(N; x, [/) is in P whenever there exists to > 1 so that 
U A {to, tv?, rnP ,...}. 

This paper is organised as previewed in this introduction. Several proofs are 
deferred to the appendix for reasons of space restriction. 

Related work. Apart from the research on circuit problems mentioned 
above there has been work on other variants like circuits over integers m and 
positive natural numbers |9], equivalence problems for circuits |16j . functions 
computed by circuits [32], and equations over sets of natural numbers [23l24j . 






2 Preliminaries 


Let N be the set of non-negative integers, and let N+ be the set of positive inte¬ 
gers. For m G N, let Div^, be the set of factors of m. Finally, let {N} be the set 
of singletons {{x} : n G N}. In this paper we use a version of the CSP permit¬ 
ting both relations and functions (and constants). Thus, a constraint language 
consists of a domain together with functions, relations and constants over that 
domain. One may thus consider a constraint language to be a first-order struc¬ 
ture. A constraint language is a core if all of its endomorphisms are embeddings 
(equivalently, if the domain is finite, automorphisms). The functional version of 
the CSP has previously been seen in, e.g., |13| . For a purely functional constraint 
language, a primitive positive sentence is the existential quantification of a con¬ 
junction of term equalities. More generally, and when relations present, we may 
have positive atoms in this conjunction. The problem CSP(T') takes as input a 
primitive positive sentence (/3, and asks whether it is true on F. We will allow 
that the functions involved on ip be defined on a larger domain than the domain 
of r. This is rather unheimlicl^ but it allows the problems of |17) to be more 
readily realised in the vicinity of CSPs. For example, one such typical domain 
is {N}, but we will allow functions such as — (complement), U (union) and fl 
(intersection) whose domain and range is the set of all subsets of N. We will also 
recall the operations of set-wise addition A + B := {a + b : a G A A b G B} and 
multiplication A x B := {ab : a G A Ab G B}. 

Sf, np, and Ap denote levels of the polynomial-time hierarchy, while Si, Bi, 
and Ai denote levels of the arithmetical hierarchy. Moreover, we use the classes 
NP = Ef, PSPACE = DSPACE(n'=), and 

3EXPSPACE = IJ DSPACE 

k>l 

For more on these complexity classes we refer the reader to |5T| . 

For sets A and B we say that A is polynomial-time many-one reducible to 
B, in symbols A<^ B, if there exists a polynomial-time computable function / 
such that for all x it holds that {x G A -4=^ /(x) G B). If / is even computable 
in logarithmic space, then A is logspace many-one reducible to B, in symbols 
Ak'^sb. a is nondeterministic polynomial-time many-one reducible to B, in 
symbols A B, if there is a nondeterministic Turing transducer M that runs 
in polynomial time such that for all x it holds that a; G A if and only if there 
exists a y computed by M on input x with y G B. The reducibility notions <Pj, 
and are transitive and NP is closed under these reducibilities. 

A circuit C = (F, E, gc) is a hnite, non-empty, directed, acyclic graph (V, E) 
with a specihed node gc G V. The graph can contain multi-edges, it does not 
have to be connected, and V = {1, 2,..., n} for some n G N. The nodes in the 
graph (F, B) are topologically ordered, i.e., for all vi,V 2 G F, if vi < V 2 , then 

Weird. Thus spake Lindemann about Hilbert’s non-constructive methods in the res¬ 
olution of Gordon’s problem (see [35]). 




there is no path from V 2 to vi. Nodes are also called gates. Nodes with indegree 
0 are called input gates and gc is called the output gate. If there is an edge from 
gate u to gate v, then we say that tt is a predecessor of v and v is the suecessor 
of u. 

Let 0 C {U, n, x}. An 0-circuit with unassigned input gates C = 
{V, E,gc, a) is a circuit {V, E,gc) whose gates are labeled by the labeling func¬ 
tion a : V ^ 0UNU{*} such that the following holds: Each gate has an indegree 
in {0,1, 2}, gates with indegree 0 have labels from N U {*}, gates with indegree 
1 have label and gates with indegree 2 have labels from {U,n,-|-, x}. Input 
gates with a label from N are called assigned (or constant) input gates; input 
gates with label * are called unassigned (or variable) input gates. An 0-formula 
is an 0 -circuit that only contains nodes with outdegree one. 

Let ui < • ■ • < Un he the unassigned inputs in C and xi,... ,Xn G N. By 
assigning value Xi to the input Ui, we obtain an 0 -eireuit C(xi,... ,Xn) whose 
input gates are all assigned. In this circuit, each gate g computes the following set 
I{g): If g is an assigned input gate where a{g) ^ *, then I{g) = { 0 ( 5 )}. Itg = Uk 
is an unassigned input gate, then I{g) = {a:^}. If g has label “ and predecessor 
gi, then I{g) = N — I{gi). It g has label o g {U,n,-|-, x} and predecessors gi 
and 52 , then I{g) = I{gi) o I{g 2 ). Finally, let I{C{xi ,... ,a:„)) = I{gc) be the 
set computed by the circuit C(xi ,..., Xn). 

Definition 1 (membership, equivalence, and satisfiability problems of 
circuits and formulas). 

Let 0C(U,n, ,-t-, x). 

MCn(0) = {{C,b) \ C is an 0-circuit without unassigned inputs and b g /(C)} 

ECn( 0) = {(Cl, ( 72 ) I Cl and C 2 are 0-eireuits without unassigned inputs and 
we have I(Ci) = /(C 2 )} 

SCn(0) = {(C, b) \ C is an 0-circuit with unassigned inputs ui < • ■ • < and 
there exist xi,... ,Xn G N such that b g /(C(a:i,... ,x„))} 

MFi!^(0), EFn(0), and SFi!^(0) are the variants that deal with 0-formulas in¬ 
stead of 0-circuits. 

When an 0-circuit is used as input for an algorithm, then we use a suitable 
encoding such that it is possible to verify in deterministic logarithmic space 
whether a given string encodes a valid circuit. 

In Section[3J for i g N, we often identify {i} with i, where this can not cause 
a harmful confusion. 


3 Circuit Satisfiability and functional CSPs 

We investigate the computational complexity of functional CSPs. In many cases 
we can translate known lower and upper bounds for membership, equivalence, 
and satisfiability problems of arithmetic circuits |28IIfilI8| to CSPs. Our main re¬ 
sult is the decidability of SCn(“, U, C, x) and CSP({N}; “, U, fl, x), which solves 




the main open question of the paper [T5] . Table [T] summarizes the results ob¬ 
tained in this section and shows open questions. In particular, we would like to 
improve the gap between the lower and upper bounds for CSP({N};0), where 
0 contains U and exactly one arithmetic operation (-1- or x). 

We start with the observation that the equivalence of arithmetic formulas 
reduces to functional CSPs. This yields several lower bounds for the CSPs. 

Proposition 1. For 0 C { , U, n, -b, x} it holds that EFn(0) CSP({N}; 0). 

Corollary 1. 

1. CSP({N}; {“, U, n,-b}) and CSP({N}; {“, U, n, x}) are -hard/or PSPACE. 
CSP({N};{U,n,+}), CSP({N};{U,n, x}), CSP({N};{U,+}), and 
CSP({N}; {U, x}) are ^^-hard for . 

CSPs with -b and x can express diophantine equations, which implies the 
Turing-hardness of such CSPs. 

Proposition 2. CSP({N};-b, x) is <J^®-hard/or Ei. 

Proposition 3. CSP({N}; U, C,-b, x) G Ei. 

Proposition 4. CSP({N}; “, U, C, -b, x) G 172. 

We show that the decidability of Skolem arithmetic m can be used to de¬ 
cide the satisfiability of arithmetic circuits without -b. This solves the main 
open question of the paper [18] and at the same time implies the decidability of 
corresponding CSPs. 

Theorem 1. SCN(",U,n, x) G 3EXPSPACE. 

Proof. Let C = C{xi,... ,Xn) be a circuit with gates gi,... ,gr, where gi,... ,gn 
are the input gates and gr is the output gate. Without loss of generality we may 
assume that C does not have H-gates. For every gate gk we define a formula (pk '■= 
(pk{xi ,..., Xn, ik,Vk, bk) in Skolem arithmetic such that the following holds. 

(*) For ai,..., a„, r G N, 6 G {0,1}, and f = 1,..., fc it holds that (pkiai, • ■ ■, On, 0, 
r, 6) is true and 

— (fikiai ,..., a„, i, V, b) is true IFF 

— (& = I iff (7(01, ■ • ■, o-n) produces at gi a set that contains v). 

Let (/?o := V -ido V {xi ■ ... ■ Xn ■ io ■ vq =0), which is always true. For 
1 < k < n, the formula ipk which corresponds to the fc-th input gate is defined 
as 

Fk '■= 34-1,Ufc_i, 6fc_i 

[(ik = k Abk = 0) ^ {xk Vk A ik-i = 0)]A 
[(ffc = fc A = I) (xfc = Vk A ik-i = 0)]A 
\fk 7^ k y (A-—1 — ik F Vk—i — Vk A hk—\ — l^fc)]A 
Fk-l- 


Observe that the free variables of tpk are variables are xi,... ,Xn,ik,Vk,bkj i.e., 

V^fc = , x„, ifc, Wfe, bk). Moreover, an induction on k shows that (*) holds 

for all (fk where 0 < k < n. 

Now define the formulas ipk for the inner gates gk where n < k < r. Here dk, 

Cfe, fk, fkj hk, and h'j, are used as auxiliary variables. 

If gk is a complement gate with predecessor gp, then let 

ifk ■= 3ik-i,Vk-i,bk-i 

[4 = fc -)> (ifc-i = p A Vk-i = Vk A (6fe = 1 bk-i = 0) A (6fe = 0 bk-i = 1))]A 
[iky^k^ (4-1 = 4 A Vk-i =Vk A bk-i = bk)] A 

If gk is a U-gate with predecessors gp and gq, then let 




3/fc, hk'iek3ik-i,Vk-i,bk-i 

[(4 = fc A e/c = 0) ->• (4-1 = p A Vk-i = Vk A 4-i = /fe)]A 

[(4 = fc A efe 0) ->■ {ik-i = q A Vk-i = Vk A 4-1 = 4)]A 

[(4 = fc A 4 = 1) (/fc = 1 V /ifc = 1)]A 

[{ik = k Abk = 0) ^ ifk = 0 Ahk = 0)]A 

[ik^k ^ (4-1 = 4 A Vk-i =Vk A bk-i = 4)]A 

‘Pfc-i- 


If gk is a x-gate with predecessors gp and pg, then let 


q^k ■ — 3fk^ fk'^^k'^hk: hf^3dk3ik—l , Vk — 1: bk— 1 

[{ik = k Abk = 1 Aek = 0) ^ {fk ■ fk = Vk A ik-i = p A Vk-i = fk A bk-i = 1)]A 

[(ifc = fc A 4 = 1 A efe 7 ^ 0) {fk ■ fk= Vk A 4-i = qA Vk-i = f'k A bk-i = ijjA 

[{ik = k A bk = 0 A hk ■ h',, = Vk A dk = 0) ^ {ik-i = p A Vk-i = hk A bk-i = 0)]A 

[{ik = k A bk = 0 A hk ■ h',, = Vk A dk ^ 0) ^ {ik-i = qA Vk-i = h'k A bk-i = 0)]A 

[ik^k ^ {ik-i = 4 A Vk-i =Vk A 4-1 = 4)]A 

Pfc-i- 

Again it holds that (pk’s free variables are xi,... ,Xn,ik,Vk,bk and an induc¬ 
tion on k shows that (*) holds for all ipk where 0 < fc < r. So for the output 
gate gr we obtain 


(C,u) € SCn( ,U, n, x)) <;=> 3ai,. 


an <Pr(ai, ■ ■ ■ ,an,r,v, I). 


The right-hand side is a first-order sentence of Skolem arithmetic. On input 
(C, v) this sentence can be computed in polynomial time, which shows that 
SCn(“, U, n, x) is <Pj-reducible to Skolem arithmetic. The latter is decidable in 
3EXPSPACE [Tg. 

Corollary 2. CSP({N}; “, U, O, x) G 3EXPSPACE 

Proof. By Theorem[Tl it suffices to show CSP({N}; “, U, 0, x) <Pj SCn(“, U, O, x). 
We describe the reduction on input of a CSP({N}; “, U, O, x)-instance x = 3j/ G 


N" Afco(^2* = t2i+i)- Observe that 


f \(Oi n t 2 i+l) U {t 2 i n t 2 i+l) — 0 

i=0 
m 

[{t2i n t2i+i) U {t2i n ^2i+l)] = 0 

i^O 

m 

0 € 0 • [(t2i n t2i+l) U (t2i n t2i+l)] ■ 

i=0 

^-V-' 

C-- 

So X S CSP({N}; U, n, x) if and only if (C, 0) € SCn(“, U, fl, x). 

Corollary 3. CSP({N}; U, n,+) G 3EXPSPACE 

Proof. By Corollary[21 it suffices to show that we have CSP({N}; U, fl, +) 
CSP({N}; U, n, x). Consider a CSP({N}; U, fl, +)-instance 

m 

x\=3y ^ N"" /\{t 2 i = Wi)- 

We may assume that 0 and 1 are the only constants that occur in x. We can 
do this, since constants c > 1 can be removed as follows: Let I = [log cj, replace 
c with a new variable z, and add constraints 

(zo = {2}) A{zi = Zq + Zq) A--- A {zi = zi^i + zi^i) A (z = ^ Zi), 

iGl 

where zq, ..., zi are new variables and 

I = {i \ the i-th bit in c’s binary representation is 1}. 

Note that removing constants in this way can be done in polynomial time. 

Observe that the term q := ^{0,1} • {0,1} fl {0,1, 2}^ • {0} fl {1} generates 
the set {2* I i G N}. 

For every term t, let t' be the term that is obtained from t if every constant 
c is replaced with 2"^, every + operation is replaced with x, and every comple¬ 
ment operation s is replaced with {s D q). The computation of t' is possible in 
polynomial time, since only the constants 0 and 1 can appear. 

The reduction outputs the CSP({N}; “, U, n, x)-instance 

m n 

x' := 3y G N” /\(4 = 4-hi) A /\( 2 /i Ug = q). 

2—0 2—1 

Observe that for all terms t and all e = (ei,..., Cn) G N" it holds that 

t'(2^C...,2'=") = {2* hGt(ei,...,e„)}. (1) 


/\ (^22 — ^ 22 + i ) 


2 = 0 













It remains to show that x and x' are equivalent. 

If e = (ei,, e„) G N" is a satisfying assignment for x, then by equation (P), 
z = (2®b ..., 2®”) is a satisfying assignment for x' (note that Ar=i(2/* U g = g) 
holds, since yi = 2®* G q). 

If z = (zi,..., z„) G N" is a satisfying assignment for x', then because of the 
constraints Ar=i(y* ^ ^ = 9); = 2®b ..., z„ = 2®” for e = (ei,..., e„) G N" 

and by (EI), e is a satisfying assignment for x. 

The following propositions transfer the NP-hardness from satisfiability prob¬ 
lems for arithmetic circuits to CSP({N}; x) and CSP({N};-f). 

Proposition 5. CSP({N}; x) is <^^-hard for NP. 

Proposition 6. CSP({N};-|-) is ^ 

The remaining results in this section show that certain functional CSPs be¬ 
long to NP. This needs non-trivial arguments of the form: If a CSP can be sat¬ 
isfied, then even with small values. These arguments are provided by the known 
results that integer programs, existential Presburger arithmetic, and existential 
Skolem arithmetic are decidable in NP |8i;i()i:i;iiiy| . 

Proposition 7. CSP({N}; “, fl, U) is complete for NF. 

Proposition 8. CSP({N}; -f) G NP. 

Proof. Consider a CSP({N}; -|-)-instance y) := 3xi,..., = fy A • ■ • A Sm = 

tm]. Each atom Si = ti term can be written as 0 = U — Si = OipXi -f • • • -I- Ui^nXn 
where Oij G Z. Let 



Hence (p G CSP({N};-|-) if and only if there exists an a: = {xi,... ,Xn) G N” 
such that Ax = 0. The right-hand side is an integer program that can be solved 

in NP [SEO]- 

Proposition 9. 


L CSP({N}; n, +) <NP CSP({N}; +, =, fy). 
CSP({N}; n, x) <NP CSP({N}; x, =, fy). 


Corollary 4. CSP({N}; C, -f), CSP({N}; C, x) G NP. 
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Table 1. Upper and lower bounds for CSP({N}; 0). All lower bounds are with respect 
to <lS®-reductions. 


4 CSPs over fo-expansions of Skolem Arithmetic 

We now commence our exploration of the complexity of CSPs generated from 
the simplest expansions of (N; x). Abandoning our set-wise definitions, we hence¬ 
forth use X to refer to the syntactic multiplication of Skolem Arithmetic (which 
may additionally carry semantic content). When we wish to refer to multipli¬ 
cation in a purely semantic way, we prefer -s or f([. We will consider x as a 
ternary relation rather than a binary function. We will never use syntactic x 
in a non-standard way, i.e. holding on a triple of integers for which it does not 
already hold in natural arithmetic. 

Proposition 10. Let F be a finite signature reduet of (N; x, 1, 2, ...). Then 
CSP{r) is in NP. 

Proof. It is known that Skolem Arithmetic admits quantiher-elimination and 
that the existential theory in this language is in NP EQ]. The result follows 
when one considers that we can substitute quantifier-free definitions for each 
among our finite set of fo-definable relations. 






















4.1 Upper bounds 

We continue with polynomial upper bounds. 

Lemma 1. Let U Q N be non-empty and U fl {0,1} = 0. CSP{N; x,U) is 
polynomial-time reducible to CSP(N~^; x,U). 

Proof. Let (p be an arbitrary instance of CSP(N; x, U) involving a set of atoms 
C on variables V. We construct a set V f-V incrementally by repeating the 
following three steps until a fixed point is reached. 

— if V{v) e U, then V' := V' U {u}, 

— if (x X y = z) € C and ^ € V, then V := V U {x, y}, and 

— if {x X y = z) G C and x,y G V, then V' := V U {z}. 

Note that if u G V', then any solution to p must satisfy s{v) yf 0. Let 
Vb = U \ V'. Construct p', an instance for CSP(N; x, U, 0), with atoms C and 
variables V' by replacing each variable v G Vq with the constant 0. Note the 
following: 

1. if 17(0) G C", then the variable that was replaced by 0 is a member of V so 
this case cannot occur. 

2. if (0 X y = z) G C' or (a; X 0 = 2 ;) G C', then the variable replaced by 0 is 
not a member of Y' while 2 is a member of Y'. This situation cannot occur. 

3. if (a; X y = 0) G C', then note that x,y gY' so the variable that was replaced 
with 0 also was a member of Y'. Hence, this case cannot occur. 

Thus, 0 can only appear in three cases: (a: x 0 = 0), (0 x x = 0), and (0x0 = 0). 
Let p" be an instance for CSP(N“'"; x, U) built from atoms C" and variables Y' 
where C" is obtained from C by removing these kinds of constraints. Note that 
(N; x,[7,0) \= p' iff (N; x,[/) |= p". Also note that if p" has a solution, then it 
has a solution s : U' —)• N'*'. This implies that p' is satisfiable on (N; x, [/) iff it 
is satisfiable on (N+; x, [/) 

The transformation above can obviously be carried out in polynomial time. In 
order to prove the lemma, it remains to show that (N; x, [/) ^ iff (N; x, 17, 0) \= 
P'- 

Assume first that p has a solution s : U N. Since C" C C, it follows 
immediately that s is a solution to p"^ too. 

Assume instead that p" has a solution s' : Y' ^ N"*". We claim that the 
function s : U —;■ N defined by s{v) = s'{v) for v G Y' and s(z;) = 0 otherwise 
is a solution to p. If the variable v G Y appears in an atom U(v), then v gY' 
and U{v) is satisfied by s. Consider an atom {x x y = z) G C.lf {x, y, z} C U', 
then s satisfies the atom since (x x y = z) G C". Assume x ^ Y'. Then z ^ U', 
s(x) = s(z) = 0 and the atom is satisfied by s. The same reasoning applies when 
y ^ Y'. Assume finally that z ^ Y'. Then at least one of x,y ^ Y' so s(x) = 0 
and/or s(y) = 0. Combining this with the fact that s(z) = 0 implies that the 
atom is satisfied by s. 


We now borrow the following slight simplification of Lemma 6 from . 

Lemma 2 (Scalability [25]). Let r be a finite signature constraint language 
with domain M, whose relations are quantifier-free definable in +, < and <, such 
that the following holds. 

— Every satisfiable instance of CSP(r) is satisfied by some rational point. 

— For eaeh relation R G E , it holds that if x := {xi,X 2 ,. ■. ,X}f) G R, then 
(axi,ax 2 ,.. ■, axk) G R for aZl a G {y : y G K., y > 1}. 

— CSP{r) is in P. 

Then CSP{A) is in P, where A is obtained from E by substituting the domain 
R by Z. 

Lemma 3. Arbitrarily choose m > 1 and U C such that ...} C 

U. Then, CSP{^'^', x,U) is in P. 

Proof. Define h(l) = 1 and h{x) = when a; > 1. Let D = {1, m, m^, m ^,... }, 
The function h is a homomorphism from (N^; x, [/) to {D; x, [/ n D). Clearly, 
h{U) = U (1 D. Suppose a-b = c where a, 6, c G N+. We see that 

h{a) ■ h{b) = = h{c). 

Define h'{l) = 0 and h'{nif) = k. Note that h' is a homomorphism from 
[D\ x,U n D) to (N; +, x > 1). We know that that CSP(R; + ,x > 0, a; > 1) is 
in P (via linear programming) and this implies tractability of CSP(N;+,x > 1) 
through CSP(Z;+,a; > l,a: > 0) by Lemma[2] 

Proposition 11. Arbitrarily choose to > 1 and U Q N such that {to, to^, to^, ...} 
C U. Then, C'6'P(N; x,U) is in P. 

Proof. Combine Lemma [1] with Lemma [S] 

4.2 Cores 

We say that an integer to > 1 has a degree-one factor p if and only if p is a prime 
such that p|to and p^ /to. Let Div^ be the set of divisors of to, pp-definable in 
(N; x,to) by 3y X X y = to. We can pp-define the relation {1} in (Div^,; x,to) 
since x = liffxxx = x (recalling 0 ^ Divm)- It follows that (1, to} are contained 
in the core of (Div^; x,m). 

Lemma 4. Let m > 1 he an integer that has a degree-one factor p. Then 
(Divm; x,to) has a two-element core. 

Proof. Consider the function e : Div^ Div^ uniquely defined by e(l) = 1, 
e(p) = TO, e(pi) = • • • = e{pk) = 1 (i.e. all the other prime divisors map to 
1), and the rule e{x ■ x') = e(x) ■ e(x'). We claim that e is an endomorphism of 
(Div(TO); x,to). Clearly, e(m) = to. Arbitrarily choose a tuple {x,y,z) G (xxy = 


z). Let X = •... • x““ and y = •... ■ y^'’ be prime factorizations. Note that 

at most one of Xi,... ,Xa,yi, ■ ■ ■ ,yb can equal p and, if so, the corresponding 
exponent must equal one. If none of the factors equal p, then e{x) = e{y) = 
e{z) = 1 and e{x) x e{y) = e{z). Otherwise, assume without loss of generality 
that xi = p. Then we have e(x) = e(z) = m and e(y) = 1. Once again e(x) x 
e(y) = e(z) and e is indeed an endomorphism of (Div^; x,m). It follows that 
({l,m}; x,m) is the core of (Div^; x,m). 

Lemma 5. Let m be an integer that does not have a degree-one factor and let 
D contain the divisors of m. Then (Div^; x,m) does not have a two-element 
core. 

Proof. Assume m has the prime factorization m = pf^ ■ ... ■ and note that 
ai,... ,ak > 1. Assume e : Div^ —> {a, b} is an endomorphism to a two-element 
core, i.e. the range of e is {l,m}. Since multiplication is determined by the 
action of the primes, we can see that for one prime p € {pi,... ,Pfe} we must 
have e(p) = m. Consider p x p = p^. If we apply the endomorphism e to this 
tuple, we end up with e(p) x e(p) = e(p)^ which is not possible. Hence, e does 
not exist and (Div^; x,m) does not admit a two-element core. 


4.3 Lower bounds 

We now move to lower bounds of NP-completeness. 

Proposition 12. CSP(N;^, x) is NP-complete. 

Proof. NP membership follows from Proposition 1101 For NP-hardness, we will 
encode the CSP of a certain Boolean constraint language, i.e. with domain {0,1}, 
with two relations: and Ri := {0,1}^\{(1,1,1)}. This CSP is NP-hard because 
^ omits constant and semilattice polymorphisms and Ri omits majority and 
minority polymorphisms (an algebraic reformulation of Schaefer’s Theorem [34] 
in the spirit of [22]). 

To encode our Boolean CSP, ensure all variables v satisfy v x v = v, which 
enforces the domain {0,1}. Consider 0 to be false and 1 to be true. 0 is pp- 
definable \>y x x x = x f\3y y ^ x f\ x x y = x. For {0,1}^ \ {(1,1,1)} take 
Ri{x,y,z) to he 3w X X y = z A w X z = 0- and for ^ take The reduction 
may now be done by local substitution and the result follows. 

An operation t : D is a. weak near-unanimity operation if t is idempo- 

tent and satisfies the equations 

t{y, x,...,x) = t{x, y,x,...,x) = ■■■ = t{x, ...,x,y). 

Theorem 2 (HI)- Let P be a constraint language over a finite set D. Lf P is 
a core and does not have a weak near-unanimity polymorphism, then CSP{r) is 
NP-hard. 


Lemma 6. Arbitrarily choose an m > I such that m ^ k'^ for all k,n > 1 
together with a finite set {l,m} C S' C N \ {0}. If (Div^; x,m) is a core, then 
CSP{S', x,m) is NP-hard. 

Proof. Assume (S; x,m) admits a weak near-unanimity operation t : ^ S. 

The relation Xi = Xk+i is pp-definable in {S; x,m) and so is the relation 

k 

R = {{xi,...,Xk) e \ = m}. 

The relation R contains the tuples 

(m, 

Applying t component-wise (i.e. vertically) to these tuples yields a tuple (a,..., a) 
for some a G D. However, R does not contain a tuple (a,..., a) for any a € D 
since this would imply that m = afi for some a,k > 1. We conclude that 
CSP(S; X, m) is NP-hard by Theorem [21 

Note that the proof of this last lemma was eased by our assumption that x is a 
relation and not a function. Had it been a function we would have to prove the 
domain S would be closed under it. 

Theorem 3. CSP{N] x,m) is NP-hard for every integer m > 1. 

Proof. If m = fc" for some k,n > 1, then we can pp-define the constant relation 
{fe} since x = k ^ nf=i ^ Hence, we assume without loss of generality 

that m k"^ for all fc, n > 1. 

We further know that Divm is pp-definable in (N; x,rn), i.e. there is polyno¬ 
mial time reduction from (Div^; x, m) to (N; x, m). The core of (Div^; x, m) is 
some (5; x, m), where {1, m} C S' C Div^ and the result follows from Lemma|6l 

Theorem 4. Let U be any subset of N\{0, 1} so that every x GU has a degree- 
one factor. Then CSPfH] x,U) is NP-hard. 

Proof. From Lemma 01 for each x G U, the core of (N; x,x) is the same (up 
to isomorphism). Fix some m GU. We claim there is a polynomial time reduc¬ 
tion from CSP(Divm; x,m) to CSP(N; x,U), whereupon the result follows from 
Theorem |3l 

To see the claim, take an instance if of CSP(Divm; x,m) and build an in¬ 
stance if of CSP(N; x,U) by adding an additional variable Vrm now substituting 
instances of m for Vm, and adding the constraint U{vm). Correctness of the 
reduction is easy to see and the result follows. 


For X G N\{0,1}, define its minimal exponent, min-exp(a;), to be the smallest 
j such that X has a factor of pi, for some prime p, but not a factor of pl^^. Thus 
an integer with a degree-one factor has minimal exponent 1. Call x S N \ {0,1} 
square-free if it omits all repeated prime factors. For a set 17 C N\ {0,1}, define 
its basis, basis(17) to be the set {min-exp(a::) : x G U}. 

Lemma 7. Let U C N \ {0,1}, so that basis(C/) is finite and basis(17) {!}. 
There is some set X pp-definable in (N; x, [/) so that basis(7f) = {!}. 

Proof. Let r = max(basis(C/)) and take an element x that witnesses this, of the 
form , where pi,... ,pk are prime and each is coprime to q (which 

is square-free), and where Oi,..., > r. Set 


f{y) := 3z,xi,... ,Xk U(z) A y"' ■ x^^ ■ ■ ■ xl’^ = z. 

We claim that f defines a set of integers X so that basis(7f) has the desired 
property. The non-emptiness is clear since 1 G basis(7f) by construction. 

Firstly, we will by contradiction argue that 0 ^ basis(7f). Assume 0 G 
basis(A). This implies that 1 G X. Hence, 3z,xi, ...,Xk U{z) AV ■ = 

z, so G U. It follows that d = min{ai,..., at} € basis(f/). Now, 

oi,..., Ofc > r so d > r. This contradicts the fact that r = max(basis(17)). 

We will now argue by contradiction that 1 < s ^ basis(Ar). Assume s G 
basis(A). Then there exists a t = q^ • "'Pk° ^ ^ where s < ai,...,ak. 
Since t G X, we know that 3z,xi,... ,Xk U{z) A P ■ = z. Let e = 

P ■ xf^ ■ • • G U as above. Let’s expand t: e = (g'’ • • • • p)!'')’’ • x“^ • • • x)!'“. 

We see that min-exp(e) > r which contradicts the choice of r. 

Example 1. We provide an example of the construction of the previous lemma in 
vivo.heiU := {p^ ,p^q^ ,p‘^q^r^ : p, g, r primes}, so that basis(f/) = {2,4}. Then 
f{y) := 3z, X U(z) A y‘^ ■ x^ = z. We can now deduce X := {p,pg^ : p, g primes} 
and basis(X) = {!}. 

Theorem 5. Let U C N\{0,1} be so thathasis{U) is finite. Then CSP{N-, x, U) 
is NP-complete. 

Proof. Membership of NP follows from Proposition [101 We use the construction 
of the previous lemma to pp-define X with basis(A) = {!}. This allows us to 
polynomially reduce CSP(N; x,X) to CSP(N; x,U) by local substitution. NP- 
hardness for the former comes from Theorem |T] and the result follows. 

5 Final remarks 

In this paper we have provided a solution to the major open question from HZ] 
as well as begun the investigation of CSPs associated with Skolem Arithmetic. 
However, the thrust of our work must be considered exploratory and there are 
two major directions in which more work is necessary. 


A perfunctory glance at the results of Section [3] shows that our bounds are 
not tight, and it would be great to see some natural CSPs in this region man¬ 
ifesting complexities such as Pspce-complete. It is informative to compare our 
Table [I] with Table 1 in m- Our weird formulation of these CSPs belies the fact 
there are more natural versions where, for O C { —, n, U, -f, x }, we ask about 
CSP(7^(N); O), where V{N) is the power set of N, rather than the somewhat 
esoteric CSP({N};0). Indeed, if we replace complement ” by set difference 
“\”, these questions could also be phrased for just the finite sets of 'P(N). These 
are all reasonable questions where we are not so sure of the boundary of decid¬ 
ability. Indeed, the complexity might be higher when the domain is 7^(N), in 
comparison to {N}. 

Meanwhile, the results of Section S] need to be extended to a classification 
of complexity for all CSP(T), where T is a reduct of Skolem Arithmetic (N; x). 
We anticipate the first stage is to complete the classification for CSP(N; x,U) 
where U is fo-definable in (N; x). We conjecture that Theorem [5] is tight in the 
sense that the outstanding cases are in P (of course this is trivial if 0 or 1 S C/). 

Conjecture 1. Let U C N\ {0,1} be fo-definable in (N; x), such that basis(C/) is 
not finite. Then CSP(N; x, [/) is in P. 

If this can be proved, the task is to extend to all problems of the form CSP(T) 
where F is an fo-expansion of (N; x). Since adding disequality already results 
in an NP-hard problem, we imagine this improvement is achievable. Finally, the 
task is to consider reducts, and not just fo-expansions, of (N; x). This might be 
harder, and will surely involve a specialised theorem of the form ‘Petrus’ in [B]. 

The mechanism through which we derived a polynomial algorithm to solve 
Proposition HD belies a close relationship between Skolem Arithmetic and Pres- 
burger Arithmetic (integers with addition). This was well-known to Mostowski 
[29] . who was able to prove decidability of Skolem Arithmetic through seeing 
it as a certain weak direct product of Presburger Arithmetic. Thus, it makes 
sense to study CSP(T), where T is a reduct of Presburger Arithmetic (N;-|-), 
in tandem with the program for Skolem Arithmetic. A first step would be to 
classify CSP(N; -f, U) where U is fo-definable in (N; -f) and a group in Dresden 
is working on this problem. Their polynomial algorithms may yet yield the 
solution to Conjectured! 
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Appendix 

Missing proofs from Section [3] 

We start with the observation that the equivalence of arithmetic formulas reduces 
to functional CSPs. This yields several lower bounds for the CSPs. 
Proposition[TJ For 0 C U, H, +, x} it holds that EFn(0) CSP({N}; 0). 

Proof. An EFN(0)-instance (^ 1 ,^ 2 ) is mapped to the CSP({N}; 0)-instance 
-El = F2. 


Corollary [TJ 

1. CSP({N}; {“, U, n, +}) and CSP({N}; {“, U, fl, x}) are <[^®-hard for PSPACE. 

2. CSP({N}; {U, n, +}), CSP({N}; {U, n, x}), CSP({N}; {U, +}), and 
CSP({N}; {U, x}) are <;;^S-hard for n^. 

Proof. The statements follow from Proposition [T] and the following facts [16] : 
EFn({~, U, n, +}) and EFN({“,U,n, x}) are <[jJ®-complete for PSPACE. 

EFn({U, n, +}), EFn({U, n, x}), EFn({U, +}), and EFn({U, x}) are <Jj]®-complete 
for pP. 


Proposition!^ CSP({N};+, x) is <|^S-hard for Ei. 

Proof. By the Matiyasevich-Robinson-Davis-Putnam theorem nna, there ex¬ 
ists an n S N and a multivariate polynomial p with integer coefScients such that 
for every A € Ai there exists an a € N such that 

xGA dy G N”,p(a, x, 1 /) = 0. 

In the equation p(a, x,y) = 0 we can move negative monoms and negative con¬ 
stants to the right-hand side. This yields multivariate polynomials I and r with 
coefficients from N such that 

xGA 3y GN'^,l{a,x,y) = r{a,x,y). 

The right-hand side is a CSP({N}; +, x)-instance. Hence A CSP({N}; +, x) 
for every A G Ai. 


Proposition [3| CSP({N}; U, fl,-b, x) G Ei. 

Proof. It is decidable whether a given assignment satisfies a CSP({N}; U, fl, +, x)- 
instance. Hence testing the existence of a satisfying assignment is in Ei. 


Proposition m CSP({N}; ,U, fl,-b, x) G E 2 . 




Proof. By Glafier et al. [TB], ECn(“, U, fl, +, x) € ^ 2 - Consider an arbitrary 
CSP({N}; U, n, +, x)-instance x := 3y € N”[to = C A • • • A {t2m = t2m+i]- It 
holds that 



The right-hand side is a IJ 2 predicate. 

Proposition [5] CSP({N}; x) is <J^S-hard for NP. 

Proof. It is known that 3SAT SCN({n, x}) [TH]. The reduction has the ad¬ 
ditional property that it outputs pairs (C, b) where the circuit C is connected in 
the sense that from each gate there exists a path to the output gate. Hence it 
suffices to construct a <|55®-reduction that works on SCN({n, x})-instances (C, b) 
where C is connected. 

For such a pair (C, b) we construct a CSP({N}; x)-instance where each gate 
g is represented by the variable g. Moreover, each gate g causes the follow¬ 
ing constraints: If g is an assigned input gate with value fc € N, then we add 
the constraint g = k. For unassigned input gates no additional constraints are 
needed. If is a x-gate with predecessors gi and g 2 , then we add the con¬ 
straint g = gi ■ g 2 - If g is a fl-gate with predecessors gi and g 2 , then we add 
the constraints g = gi and g = g 2 . li g is the output gate, then this causes 
the additional constraint g = b. Finally, if gi,...,gn are the gates in C and 
Cl,..., Cm are the constraints described above, then the reduction outputs the 
CSP({N}; x)-instance ip := 3gi, ..., g„[ci A • • • A Cm]- 

It remains to argue that for connected C it holds that 


(C',&)eSCN({n,x}) ^ ^gCSP({N};x). 


Assume (C, 6) G SCN({n, x}) and consider an assignment that produces {b} 
at the output gate. Since C is connected, each gate gi computes a singleton { 0 ^}. 
Hence oi,..., a„ is a satisfying assignment for ip, which shows ip G CSP({N}; x). 

Assume ip G CSP({N}; x). Let oi,... ,a„ be a satisfying assignment for ip 
and let I be the number of C’s input gates. The constraints in ip make sure that 
(7(01, ■ ■ ■ ,0,1) produces {oi} at gate gi. In particular, (7(ai, ... ,ai) produces {6} 
at the output gate, which shows {C, b) G SCN({n, x}). 

Proposition [6] CSP({N}; -I-) is <(((®-hard for NP. 

Proof. It suffices to show SCn({-|-}) <m®CSP({N}; -b) [H]. The proof is similar 
to the proof of Proposition [SJ but easier, since we have no fl-gates and hence we 
do not need the assumption that C is connected. 


Proposition [3 CSP({N}; , H, U) is <(((S-complete for NP. 


Proof. Consider a CSP({N}; Cl, U)-instance ip := 3xi,..., A • • • A 

tm = im]- We show that if € CSP({N}; Cl, U), then it has a satisfying 
assignment a = (oi,..., a„) such that ai,..., G {0,..., n — 1}. This implies 
CSP({N};-,n,U) G NP. 

Assume ip G CSP({N}; Cl, U) and choose a satisfying assignment a = 
(oi ,... ,an) such that a' = max{ai,..., a„} is minimal. Assume that a' > n, 
we will show a contradiction. Let b' = min(N — {ai,...,a„}) and note that 
b' < n. Let b be the assignment that is obtained from a if all occurrences of 
a' are replaced with b'. For any term t in (p, the sets computed by t under the 
assignments a and b are denoted by fy and fy, respectively. Observe that for all 
terms t in (p and all a; G N — {o', 6'} it holds that: 

X € ta X € tb 
q! G ta ) -• b' G tb 
y G ta ) -• q! G tb 

It follows that for all atoms f = f' in it holds that 


Therefore, & is a satisfying assignment that is smaller than a, which contradicts 
the minimal choice of a. 

For the NP-hardness it suffices to show 3SAT CSP({N}; fl, U). On 
input of a 3CNF-formula t = t(xi,... ,Xn) the reduction outputs the instance of 
CSP({N};-,n,U) 

ip := 3xi,.. .,x„[t' n {1} = {!}], 

where t' is obtained from t by replacing -i, A, V with O, U, respectively. Every 
satisfying assignment for t also satisfies p. Conversely, if ai,..., On is a satisfying 
assignment for (p, then we obtain a satisfying assignment for t if values greater 
than 1 are replaced with 0. 


Proposition [9] 

1. CSP({N}; n, +) <NP CSP({N}; +, =, fy). 

2. CSP({N}; n, x) <NP CSP({N}; x, =, fy). 

Proof. We show the first statement, the proof of the second one is analogous. 

For a term t, let t' be the term obtained from t if every subterm of the form 
Si n S 2 is replaced with si. 

We describe the <^^-reduction on input of a CSP({N}; fl, +)-instance 

P .— 3^:1, . . . , Xyi [to — 11 A ■ • ■ A t2m — t2m+l\ • 

For each atom t 2 i = t 2 i+i, we guess nondeterministically whether t 2 i = t 2 i+i G 
{N} or t 2 i = t 2 i+i = 0. If we guessed t 2 i = t 2 i+i G {N}, then replace t 2 i with 
t 2 i, replace t 2 i+i with for every subterm si fl S 2 that appears in t 2 i 


or t 2 i+i add the constraint = s^. If we guessed ^ 2 * = t 2 i+i = 0, then guess 
a subterm ui fl U 2 in t 2 i, guess a subterm U 3 n U 4 in t 2 i+i, remove the atom 
t 2 i = i 2 i+i, and add the constraints u[ ^ U 2 and ^ U 4 . The obtained formula 
'(/' is the result of the <^^-reduction. 

We argue that the described -reduction reduces the CSP({N}; n,-b) to 
theCSP({N};+,=,^). 

Assume (/? G CSP({N}; n, -b) and fix some satisfying assignment a = (oi,..., 
an) & N". Consider the nondeterministic path of the reduction that for all atoms 
correctly guesses whether t 2 i = t 2 i+i G {N} or t 2 i = t 2 i+i = 0, and that for all 
t 2 i = t 2 i+i = 0 guesses subterms ui n U 2 in t 2 i and M 3 fl U 4 in t 2 z+i such that 
ui,U 2 , M 3 , M 4 G {N}, Ml ^ U 2 , and M 3 ^ U 4 . If t 2 i = t 2 i+i G {N}, then = hi = 
hi+i = ^ 21+1 hence the formula is still satisfied after replacing t 2 i with 
and t 2 i-i-i with Moreover, the added constraints = S 2 are satisfied, since 

in t 2 i and hi+i all subterms si fl S 2 must be nonempty. If hi = hi+i = 0, then 
after removing the atom t 2 i = hi+i and after adding the constraints u[ ^ m^ 
and M 3 7 ^ M 4 the formula is still satisfied. So at the described nondeterministic 
path the reduction outputs a formula ^ G CSP({N}; -b, =, 7 ^). 

Assume there is a nondeterministic path where the reduction outputs a for¬ 
mula Ip G CSP({N}; -b, =, 7 ^). Consider a satisfying assignment a for ip, we claim 
that a satisfies tp. If this is not true, then ip must have an atom t 2 i = hi+i that 
is not satisfied by a. 

Case 1: At the path that produced ip we guessed that hi = hi+i G {N}. 
In this case we added the constraints which ensure that hi = ^21 

hi+i = t 2 i+i- Hence under the assignment a it holds that hi = ^21 = ^ 2 i+i = 
hi+i, which contradicts the assumption that hi = hi+i is not satisfied by a. 

Case 2: At the path that produced ip we guessed that hi = hi+i = 0- Here we 
added the constraints u'^ 7 ^ u '2 and M 3 7 ^ M 4 , which are satisfied by a. Hence under 
the assignment a we have hi = hi+i = 0 , which contradicts the assumption that 
hi = hi+i is not satisfied by a. 

It follows that a satisfies p and hence p G CSP({N}; C, -b). 


Corollary H CSP({N}; fl,-b), CSP({N}; C, x) G NP. 

Proof. CSP({N}; -b, =, 7 ^)-instances and CSP({N}; x, =, 7 b)-instances are formu¬ 
las of existential Presburger arithmetic and existential Skolem arithmetic, which 
are both decidable in NP |33ll9j . Now the statement follows from Proposition [HI 




